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ABSTRACT 



Aims. X-ray properties of the stellar population in the Carina OB 1 association are examined with special emphasis on early-type 

stars. Their spectral characteristics provide some clues to understanding the nature of X-ray formation mechanisms in the winds of 

single and binary early-type stars. 

Methods. A timing and spectral analysis of five observations with XMM-Newton is performed using various statistical tests and 

thermal spectral models. 

Results. 235 point sources have been detected within the field of view. Several of these sources are probably pre-main sequence stars 

with characteristic short-term variability. Seven sources are possible background AGNs. Spectral analysis of twenty three sources of 

type OB and WR 25 was performed. We derived spectral parameters of the sources and their fluxes in three energy bands. Estimating 

the interstellar absorption for every source and the distance to the nebula, we derived X-ray luminosities of these stars and compared 

them to their bolometric luminosities. We discuss possible reasons for the fact that, on average, the observed X-ray properties of 

binary and single early type stars are not very different, and give several possible explanations. 
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1 . Introduction 

Although X-ray emission by early-type stars is nowadays well 
established, its origin is still not fully understood. Observations 
with the Einstein, i^OSAT and ASCA satellites indicated that hot 
luminous stars have rather soft thermal X-ray spectra. This pic- 
ture contrasts with the expected spectral properties if the X-rays 
were produced in a hot corona at the base of the stellar wind (e.g. 
Waldron il98 4). Indeed, in the latter case one would expect to ob- 
serve substantial absorption of the softest emission. The lack of 
this absorption in the observed spectra triggered the elaboration 
of an alternative scenario. In fact, according to the phenomeno- 
logical model proposed by Lucy & White d 19801 1 and further 
elaborated by Lucy ( j 19821 ). hydrodynamic shocks are generated 
throughout a radiation driven stellar wind as the consequence of 
the intrinsic instabilities of radiative driving. The velocity jumps 
in such shocks heat the post-shock plasma to temperatures of a 
few million degrees. These shocks would be distributed through- 
out the stellar wind with the result that even soft X-rays could es- 
cape the wind without substantial absorption. Recent theoretical 
studies of such instabilities have been presented for instance by 
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Feldmeieretal. ( fT997b . Owocki & Cohen (T9W), and Dessart & 
Owocki ( !2002l l. High resolution spectra obtained with the new 
generation of X-ray observatories have confirmed the thermal 
origin of the bulk of the X-ray emission from early-type stars 
(e.g. Kahn et al. |2001| l. but in some cases, the properties of the 
observed line profiles led also to new challenges for the shock 
model (see e.g. Rauw ' 20061 for a review). 

Another well established, though yet unexplained, property 
of the X-ray emission of hot stars is the linear scaling be- 
tween the observed X-ray luminosity (Lx) and the bolometric 
luminosity (Lboi)- Using Einstein data. Long & White (1980|, 
Pallavicini et al. ( 11981b and Chlebowski et al. ( 1989) found that 
Lx — 10~^ X Lboi for OB stars. ROSAT observations broadly 
confirmed this relation. Berghofer et al. ( |19971 l and Kudritzki et 
al. ( 1996) found that including a weak dependence on the char- 
acteristic wind density leads to a somewhat tighter relationship. 
However, considering the scaling of the X-ray emission with var- 
ious wind and stellar properties in the framework of the shock 
model, Owocki & Cohen ( 119991 ) found that a delicate equilib- 
rium between emission and absorption is required to reproduce 
the empirically derived Lx - ^boi relation. With the broader en- 
ergy range and improved sensitivity of the new generation of 
X-ray satellites it is important to re-address the question of the 
empirical Lx - iboi scaling for OB stars. 

When doing this it is crucial to pay attention to the multiplic- 
ity of the stars in the sample. Indeed, in OB binaries, the collision 
of the two stellar winds is expected to generate a strong X-ray 
bright shock between the stars (see e.g. Stevens et al. 119921 ) and 
this feature is believed to manifest itself as an "excess" in the 
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Lx/iboi ratio (see e.g. Pollock fl 9951 ). Moreover, since the rel- 
ative velocities of colliding winds can be much higher than the 
shock velocity jumps in winds of single stars, the X-ray spectra 
of colliding wind early type binaries can be quite hard. In fact the 
hardness of the observed spectra of some early-type stars, and, 
in particular, the presence of the Fe XXV and Fe XXVI lines at 
6.7keV indicative of very hot material has been suggested as an 
indication of possible binarity (e.g. Raassen et al., 2003). 

This latter suggestion leads to an interesting question. 
Clearly in a general case it is very difficult to "hide" such a hard 
secondary component of a colliding wind binary. Indeed, to pro- 
duce strong hard X-ray emission, the companion must be a mas- 
sive hot star, and as such, should manifest itself in other wave- 
length domains (e.g. optical spectra) and also via orbital spectral 
and photometric variability in the X-rays. Traditionally, a lack 
of variations was attributed to a "pole-on" orientation of the bi- 
nary orbit. However, lately, a growing number of examples of 
apparently single O and WR stars with hard X-ray spectra (e.g. 
Skinner et al. 120021 Raassen et al. 120031 the present study) has 
become available. It seems unlikely that all these objects are bi- 
naries seen from their orbital poles. Thus the question is whether 
there exists another, yet unknown, mechanism which could pro- 
duce hard X-rays in single star winds? Exceptions to this rule 
are stars such as 9^ Ori C with strong enough magnetic fields to 
confine the stellar wind into the plane near the magnetic equator 
(Cassinelli et al., 2003, Bro wn et a l., 2003, Maheswaran 2003, 
Gagne et al. 2005 Li et al., 120071 ). In such cases, the head-on 
collision of the winds from the two hemispheres of a single star 
produces a rather hard X-ray emission. The fundamental prob- 
lem here is that the only source of energy in the classical pic- 
ture of a hot star wind which is able to produce high tempera- 
tures and hence hard thermal X-rays, is the wind kinetic energy. 
But to release most of this energy through radiation one needs 
a "wall" able to stop the supersonic flow and heat the material 
in the formed shock. In the winds of single stars we do not ex- 
pect to find such a wall except in the magnetised case. On the 
other hand, there is a number of known early-type binaries which 
show neither significant X-ray excess nor particularly hard X-ray 
spectra. What could be the mechanism which dumps the X-ray 
production in these systems? 

The Carina OB 1 association with its impressive number of 
early-type stars - including some of the youngest, hottest and 
most massive objects known in our Galaxy - is an ideal place to 
investigate the X-ray properties of these objects. The history of 
X-ray studies of the region goes back to the Einstein era. Seward 
et al. (119791 ) presented Einstein X-ray (0.2-4.0keV) observa- 
tions of the Carina open cluster Tr 16 and its environment, in- 
cluding six O-type stars and one WR star, WR25. Subsequent 
Einstein observations by Seward & Chlebowski (1982) of the 
same region showed X-rays from 15 O-type and WR stars. More 
studies were performed with the use of the ROSAT and ASCA 
satellites. Some of them were focused on individual objects like 
the LBV ryCar (Corcoran et al. 2000 and references therein) or 
WR25 (e.g. Pollock & Corcoran 2006, Skinner et al. [T995] ). 
Others were devoted to investigating general X-ray properties of 
the early-type stars population (Corcoran 119991 and references 
therein). 

Recently, four studies of the region based on X-ray obser- 
vations with Chandra and XMM-Newton satelUtes were pub- 
lished. Evans et al. (120031) (hereafter EV03) and Evans et al. 
020041 ) (hereafter £^04) using Chandra data, detected 154 point 
sources, 23 of which are of O, B types. They present luminosi- 
ties and hardness ratios for the detected sources and confirm the 
Lx - iboi relation for early-type stars. They also discuss the 



low-resolution spectra of the 14 brightest sources. Sanchawala 
et al. ( 120071 ), using the same Chandra data as EV03 plus an- 
other archival Chandra data set, detected 454 sources in the re- 
gion, among which 38 are known OB stars. The latter paper is 
focused on cross-identification of the X-ray sources with opti- 
cal and near-infrared imaging observations specifically obtained 
for this study. This allowed the authors to conclude that about 
300 sources are likely late-type pre-main-sequence stars. The X- 
ray luminosities of the detected OB stars were estimated based 
on their count rate and using a single-temperature Raymond- 
Smith model of thermal plasma with log T = 6.65 (0.384 keV). 
Albacete-Colombo et al. (2003) (hereafter AC03) detected 80 
discrete sources in the region using two observations by XMM- 
Newton (revolution numbers 115 and 116). They discuss X-ray 
properties of the sources, including the Lx - iboi relation for 
early-type stars. In the present paper, we used three additional 
XMM-Newton observations, which allow us to increase the sig- 
nal to noise ratio and to detect more than 200 discrete sources. 
We shall compare our results with those of AC03 throughout 
the current paper. We do not present a detailed spectral analy- 
sis of individual objects in this paper as many such objects de- 
serve dedicated papers. A detailed study of WR25 based on the 
same XMM-Newton data as in the current study was published 
by Raassen et al. (120031 ). whilst 77 Car has been the subject of 
many studies including those utilizing extensive data sets ob- 
tained with various X-ray observatories (see, e.g., Leutenegger 
et al. l2003t Weis etal.l2504: Viotti et al. l2004l Corcoran 2002 
Hamaguchi et al. 120071 ). Analysis of other interesting objects 
(e.g. the early type binary HD 93205) will be the subject of forth- 
coming papers. 

The goals of this paper are: 

- To describe general X-ray properties of the detected sources. 

- To look for their variability, both short and long term. 

- To investigate the Lx/^boi relation for early type stars. 

- To compare X-ray properties of single versus binary early- 
type stars and to try to get some clues to their differences 
and similarities. 

The structure of the paper follows these outlined priori- 
ties. In the second section, the observational details are given. 
Source detection, their identification with the optical and IR 
catalogues. X-ray count rates, identification of possible extra- 
galactic sources are given in Section 3. Section 4 provides gen- 
eral information about open clusters in the Carina nebula, along 
with the discussion of their distance from the Earth. The results 
of the variability study are presented in Section 5. The X-ray 
properties of early-type stars detected in the region are presented 
in Section 6. We discuss our results in section 7 and a short sum- 
mary is given in Section 8. 



2. Observations 

The observations of the Carina region were obtained with the 
Reflection Grating Spectrometers (RGS) and the European Photon 
Imaging Camera (EPIC) CCD detectors of the XMM-Newton ob- 
servatory. The log of EPIC observations is shown in Table [T] In 
the two first data sets the primary target was 77 Car while in the 
third to fifth data sets the primary target was WR25. The anal- 
ysis of RGS and EPIC observations of WR 25 was presented in 
Raassen et al. (2003) while the RGS analysis of 77 Car was pub- 
lished by Leutenegger et al. (2003) and we will not repeat these 
here. In Fig. [T] we show the mozaiced image of the area com- 
bined from all 5 data sets and 3 EPIC instruments. 
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Table 1. Log of the Carina region observations by XMM- 
Newton. 



Data Set Number 
Revolution 
Obs. Date 


12 3 4 5 

#115 #116 #283 #284 #285 

26-7-2000 27-7-2000 25-6-2001 28-6-2001 30-6-2001 


Start (UT) 


04:58 


23:48 


06:51 


07:22 


04:38 


Instrument 




Inte 


gration Time 
(hr) 




MOS 1 
MOS 2 
PN 


9.4 
8.5 
8.8 


3.1 
2.3 
2.6 


10.2 
10.2 
9.6 


11.7 
11.7 
11.0 


10.4 
10.4 
9.7 



The EPic-MOS and pn instruments were operated in the full 
frame mode except during rev. 115 and 116 where MOS2 was 
operated in the small window mode. All three EPIC instruments 
used the thick filter to reject optical light. We used version 7.0 
of the XMM-Newton Science Analysis System (SAS) software 
to reduce the raw EPIC data. For spectral analyses, we adopted 
the most up-to-date redistribution matrices provided by the EPIC 
instrument teams and used SAS to build the appropriate ancillary 
response file for each observation. More details on the pipeline 
processing of the data are given in Raassen et al. (I2003I I. 

3. Source detection and cross-identification 

3.1. Source detection 

Sources were detected and count rates measured in three energy 
bands - soft (0.4- l.OkeV), medium (1.0-2.5 keV), and hard 
(2.5- lO.OkeV). The choice of these bands is based on the fol- 
lowing considerations: 

- There is a lot of noise (sometimes systematic features) in the 
PN data fori? < 0.4 keV. 

- The EPIC sensitivity above lOkeV is almost zero. 

- The relatively narrow soft band is sensitive to interstellar ab- 
sorption. 

- The wide hard band allows one to collect more photons and 
hence improve statistics. 

Source detection and determination of source parameters 
was performed with the SAS edetect.chain metatask based 
on the sliding cell detection and maximum likelihood (ML) 
methods. The ML method (the SAS task emldetect) makes 
use of maximum likelihood PSF (point spread function) fitting 
to the source count distribution. The so-called logarithmic de- 
tection likelihoods obey the simple relationship L2 — — ln(p) 
where p is the probability for a random Poissonian fluctuation to 
have caused the observed source counts. The threshold for inclu- 
sion of sources in the final output was set to i™'" = 10.0. In the 
limit of a very large number of counts, the likelihood function 
is a gaussian and the logarithmic likelihood becomes a parabola. 
At this limit the value of 10.0 is equivalent to 3 sigmas. 

To increase the signal-to-noise ratio and thus to detect as 
many sources as possible, we merged the event lists and images 
from different data sets. Unfortunately, the angular distance be- 
tween the central axes of the data sets 1 and 2 on one hand, and 
3 to 5 on the other, is too large (^ 7') and does not allow one to 
merge all data sets into a single coordinate system without pro- 
ducing spurious results in the subsequent source detection. For 



this reason we merged the data sets 1 and 2, and the data sets 3, 4, 
5 separately. The source detection was performed on the two re- 
sulting merged sets of data, making use of all energy bands and 
all instruments simultaneously. The resulting source lists were 
visually inspected and a few spurious detections, mainly along 
the CCD edges, were removed. 

The coordinates of the sources common in the two merged 
data sets are very close and do not show any systematic differ- 
ences. We attempted to use the eposcorr procedure to im- 
prove the coordinates by cross-correlating the XMM coordinates 
with the Two Micron All Sky Survey point source catalogueQ 
(2MASS, Skrutskie et al. 1997) and the optical sources in the 
SIMBAD database. In both cases, neither a significant system- 
atic shift nor a field rotation was apparent. The corrections in RA 
and DEC were less than 0.5"and not consistent between the IR 
and optical catalogues. We conclude that the procedure does not 
allow one to improve the XMM-Newton coordinates. For the fi- 
nal source list the coordinates of the sources common in the two 
merged data sets were averaged. As some of the detected sources 
may be variable, it is of interest to measure their count rates not 
just in the merged data sets but in every individual one. For this 
purpose, we repeated the ML routine with our source list as the 
input, for the data sets 1-5. 

The total number of detected sources is 235. Recall that 
AC03, using the data sets 1 and 2, detected 80 sources. 
Sanchawala et al. (2007) were able to detect 454 sources in the 
Chandra data on the field, thanks to the superior angular resolu- 
tion of the Chandra telescope. FigHI shows the positions of the 
detected sources over the mozaiced image of the field. A sample 
of the catalogue of the detected discrete X-ray sources is pre- 
sented in Table|2] The full table is available as online data. The 
table includes only those measurements in which the total loga- 
rithmic EPIC likelihood exceeds the value of 10.0. 

One should note that the count rates in Table|2]may be sys- 
tematically different from those provided for the same stars by 
AC03. One evident reason is that we used different energy lim- 
its for our soft, medium, and hard bands than AC03. A deeper 
reason, is, however, that AC03 used the wavelet technique (the 
SAS ewavelet routine) to estimate their count rates. Using 
ewavelet is equivalent to fit a Gaussian to the observed dis- 
crete source image. However, the actual EPIC PSF is quite differ- 
ent from a Gaussian. The emldetect routine makes use of an 
empirical PSF and thus, according to the SAS manual, should 
provide more accurate count rate estimates than ewavelet. 
Our tests show that the AC03 count rates are systematically 
lower than the count rates obtained for the same energy band 
with emldetect. 

Fluxes and luminosities can be calculated from count rates 
provided that the response function of the telescope is known 
and certain assumptions are made about the emission model, in- 
terstellar absorption and distance to the region. We used xspec 
to calculate the conversion factors. A combination of models 
wabs*apec was used. Here apec is an emission spectrum 
from a hot optically thin thermal plasma. The APEC code (Smith 
et al., I2OOII 1 is a modern version of the well-known model 
of Raymond & Smith (19771). wabs is an interstellar medium 
(ISM) absorption model. As we are most interested in X-ray 
properties of hot massive stars, the plasma temperature in the 
apec model was set to kT = 0.6 keV, a value typical for X-ray 
emission produced by shocks formed in stellar winds of early- 
type stars. The interstellar absorption over the field is relatively 
uniform with £;(i3-V^) = 0.52 mag (Massey & Johnson. il993l l. 
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Table 2. Sample of the X-ray catalogue. 



X# Var Status Data Set 

[1] [2] [3] 



Uncert 

Const 

Noinfo 
Var 

Uncert 

Const 
Const 
Const 



X# Var Status 

[1] [2] 



Data Set 

[3] 



1 Uncert 

2 Const 

3 Noinfo 

4 Var 



Uncert 

Const 
Const 
Const 



X# Var Status Data Set 

[1] [2] [3] 



1 Uncert 

2 Const 

3 Noinfo 

4 Var 



Uncert 

Const 
Const 
Const 



rPN 

1^2 CrpN 



(JpN 



[4] 



[5] 



[6] 



cr-pN 



[7] 



"~s — 

CTpN 



[8] 



crpN 



[9] 



CTpN 

[10] 



CpN 

[11] 



—R — 

(TpN 

[12] 



37.0 

67.8 

78.9 

155.1 

5.9 

14.2 

257.1 

165.1 

27.6 

46.6 

25.8 

41.6 

59.4 

27.9 

3.0 

64.3 

81.9 

89.4 



11.6964 

14.3821 

11.0794 

13.0312 

1.5988 

13.1136 

30.2308 

18.3160 

6.2785 

8.9067 

4.7296 

8.5184 

8.3989 

6.7066 

2.7511 

11.2433 

11.4430 

11.0901 



1.8388 
1.6988 
1.5107 
1.2820 
0.6920 
3.7333 
2.6712 
1.7907 
1.2448 
1.4053 
0.9220 
1.6391 
1.3907 
1.7137 
1.4317 
1.4168 
1.3621 
1.2253 



7.5714 
9.8716 
1.0840 
0.0000 
0.0000 
3.0763 
6.3221 
5.1892 
4.3118 
5.4440 
3.1683 
0.7818 
0.0000 
0.9344 
0.7435 
7.0934 
6.0041 
7.4002 



1.5446 
1.4067 
0.7725 
0.3180 
0.3058 
2.6452 
1.5744 
1.2056 
0.9343 
1.0365 
0.7315 
0.9702 
0.7206 
1.1214 
0.9801 
1.1558 
1.0333 
1.0139 



4.1250 
3.6847 
4.9085 
5.0587 
0.0631 
6.8288 
19.5985 
10.6765 
1.2414 
2.8946 
1.3936 
2.2362 
2.0239 
5.7722 
2.0076 
4.1499 
5.4389 
3.6899 



0.9757 
0.7888 
0.8522 
0.7272 
0.2185 
2.0792 
1.8292 
1.1067 
0.5292 
0.7262 
0.4500 
0.7546 
0.6175 
1.1112 
0.8730 
0.7856 
0.8331 
0.6185 



0.0000 
0.8258 
5.0868 
7.9724 
1.5357 
3.2085 
4.3102 
2.4504 
0.7254 
0.5681 
0.1676 
5.5004 
6.3750 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 



riviuyi 

L2 crtAOSi 



CTMOSl 



[13] 



[14] 



[15] 



CmOSI 



[16] 



"~s 



[17] 



C^MOSl 



[18] 



M- 

''"MOSl 



[19] 



[20] 



2.8 

30.0 

41.7 

54.6 

6.3 

26.8 

116.8 

66.2 

10.7 

25.3 

6.5 

8.0 

15.4 

3.8 

1.5 

26.0 

39.5 

27.0 



1.5736 
4.5064 
3.7951 
4.0763 
1.1639 
5.4286 
11.9558 
7.3132 
1.8870 
3.5911 
1.2127 
1.8536 
2.5796 
1.6451 
0.9795 
3.5804 
4.2651 
3.1810 



0.7445 
0.8380 
0.7415 
0.6695 
0.4805 
1.0912 
1.3586 
1.0267 
0.6134 
0.8622 
0.4765 
0.7433 
0.7105 
0.7503 
0.6168 
0.7153 
0.7098 
0.5963 



0.2249 
2.3016 
0.0000 
0.1027 
0.0000 
1.3413 
2.6415 
1.9323 
0.7426 
1.1051 
1.0220 
0.0000 
0.1040 
0.2117 
0.0000 
1.7610 
2.9362 
1.9396 



0.3832 
0.6342 
0.1756 
0.2009 
0.1250 
0.6648 
0.7379 
0.6414 
0.3567 
0.4992 
0.4014 
0.3138 
0.2827 
0.4138 
0.2356 
0.5240 
0.5725 
0.4694 



1.1755 
2.2048 
2.1814 
1.7207 
0.4476 
3.0680 
6.3675 
3.8654 
0.8337 
2.4859 
0.1907 
0.2292 
0.8577 
0.9603 
0.6828 
1.7545 
1.0653 
1.2414 



0.5769 
0.5365 
0.5132 
0.3944 
0.3061 
0.7086 
0.8938 
0.6433 
0.3534 
0.6274 
0.2387 
0.3111 
0.4188 
0.4802 
0.4627 
0.4686 
0.3545 
0.3615 



0.1732 
0.0000 
1.6137 
2.2529 
0.7163 
1.0192 
2.9468 
1.5155 
0.3107 
0.0000 
0.0000 
1.6244 
1.6178 
0.4731 
0.2968 
0.0648 
0.2635 
0.0000 



"TTHUHT" 

[22] 



[23] 



[24] 



[25] 



[26] 



[27] 



[28] 



[29] 



16.9 
21.8 
16.9 
19.2 

5.3 
35.8 
75.6 
45.2 

1.9 
34.9 
14.5 

7.7 
15.7 

3.6 

1.7 
22.2 
22.4 
18.6 



3.7830 
3.3101 
3.7571 
3.6836 
1.2067 
5.1319 
7.5895 
4.1332 
0.6367 
3.3469 
1.7518 
1.3382 
1.9042 
1.7374 
0.6576 
3.4013 
2.9989 
2.6419 



0.9184 
0.7358 
1.0486 
0.9986 
0.5056 
0.9696 
1.0367 
0.7263 
0.5123 
0.6648 
0.4607 
0.6299 
0.5742 
0.7574 
0.5783 
0.7177 
0.6168 
0.5682 



2.4310 
1.2515 
0.0000 
0.0000 
0.0287 
1.5404 
1.4724 
0.2420 
0.6367 
1.6299 
1.0534 
0.0000 
0.0000 
0.6287 
0.0000 
2.2301 
1.7462 
1.4887 



0.7169 
0.5254 
0.2600 
0.1867 
0.1578 
0.6497 
0.5831 
0.3470 
0.3823 
0.4648 
0.3500 
0.3383 
0.2060 
0.5433 
0.3803 
0.5873 
0.4723 
0.4440 



1.3520 
1.9346 
1.8097 
0.8314 
0.3144 
2.7681 
4.3611 
2.9420 
0.0000 
1.6700 
0.6985 
0.0900 
0.4021 
0.9286 
0.6576 
1.1712 
1.2527 
1.1532 



0.5082 
0.4799 
0.6987 
0.5721 
0.2606 
0.6045 
0.7020 
0.5314 
0.1562 
0.4481 
0.2908 
0.2523 
0.3204 
0.4724 
0.3924 
0.4071 
0.3939 
0.3414 



0.0000 
0.1240 
1.9474 
2.8522 
0.8636 
0.8234 
1.7560 
0.9493 
0.0000 
0.0470 
0.0000 
1.2483 
1.5021 
0.1802 
0.0000 
0.0000 
0.0000 
0.0000 



0.2079 
0.5338 
0.9793 
1.0067 
0.5810 
1.6178 
1.1448 
0.7269 
0.6297 
0.6108 
0.3355 
1.0844 
1.0165 
0.6666 
0.5718 
0.2330 
0.3058 
0.3015 



C^MOSl ''"MOSI 



[21] 



0.2732 
0.1103 
0.5057 
0.5023 
0.3486 
0.4967 
0.7088 
0.4785 
0.3524 
0.3171 
0.0948 
0.5977 
0.4994 
0.4013 
0.3328 
0.1322 
0.2246 
0.0678 



CrMOS2 CrMOS2 CrMOS2 <'"mOS2 CrMOS2 ''"MOS2 Cr]y[oS2 ''"MOS2 



[30] 



0.2671 
0.1872 
0.7375 
0.7969 
0.4035 
0.3904 
0.4919 
0.3531 
0.3032 
0.1584 
0.0721 
0.4676 
0.4297 
0.2353 
0.1893 
0.0666 
0.0468 
0.0958 



The first column gives the source number. The second column provides the set-to-set variability status of the source (see section 5.2). The third 
column gives the data set number in which the count rates of the source were measured. Cols. 4-12 (resp. 13-21 and 22-30) give the logarithmic 
likelihood L2 for the given instrument, the total count rate cr in the whole energy band (0.4-10.0 keV) and its associated error (a), the count rates 
in the three different energy bands (S : [0.4 - 1.0 keV], M : [1.0 - 2.5 keV], H : [2.5 - 10.0 keV]) and their errors. The count rates and the 
related uncertainties are all expressed in 10~^ counts s~ . 
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Fig. 1. Mosaiced exposure-corrected image of the field, combining all data sets and all instruments MOS 1 , MOS2, and PN. Energy 
range 0.4 - lO.OkeV. Positions of the detected sources are shown. 



Using a well-known transformation from Bohlin et al. d 19781 1 
Nh = 5.8 X 10^^ xE{B — V) this translates to a column density 



Nh 



10 cm . We used this value in the wabs model. 



Using the above combination of models and XMM-Newton re- 
sponse function we derive the conversion factors for MOS and 
PN instruments in the 0.4- lO.OkeV energy range as 

Fa, MOS = 0.76 X 10""ergscm"^s"^ x cr 

Fa,pN = 0.23 X 10""ergscm"^s"^ x cr 

Fn MOS = 2.01 X lO^^erffscm^^s"^ x cr 



O.MOS 



0,PN 



2„-l 



0.61x10 ergs cm "s x cr 



where cr is the corresponding count rate, the indices "a" and 
"0" mark the absorbed and unabsorbed fluxes respectively. 
Luminosities can be readily obtained from the above fluxes as- 
suming a distance to the nebula of 2.5 kpc (see section 6). We 
caution that the conversion from count rates to energy is very 
model-dependent and can only be used as a very rough esti- 
mation of the flux. Below we provide individual fluxes for the 
brightest sources which allow spectral fitting. 

3.2. The detection limit 

Evaluation of the detection limit of our observations is not a 
trivial task. First, this limit is a priori not uniform throughout 
the field of view as the XMM-Newton effective exposure dura- 
tion is decreasing from the center of the FOV towards the edges. 



Second, as different sources in our catalogue were detected in 
different observations or their combinations, the detection limits 
in the intermediate partial source lists are different. 

Thus, to estimate the overall limit, we adopted a completely 
empirical approach, taking advantage of the large number of X- 
ray sources in the field. We assumed that a good indication of 
the detection limit in the different parts of the field is given by 
the brightness of the faintest sources detected in these areas. Due 
to the presence of gaps, not all sources which are detected in the 
MOS images, are also present in the PN ones. To increase statis- 
tics, we first computed an equivalent PN count rate, in the range 
0.4-10. OkeV, for each source missing in the PN images. To first 
order, the relation between the count rates measured in any of 
the two MOS detectors and in the PN detector is linear. We found 
the conversion factor using the count rates of the sources, which 
were detected in all three EPIC instruments. Figure|2]displays the 
source PN-equivalent count rates as a function of the distance 
from the central axis of the FOV. A lower limit is clearly seen 
in the distribution. Towards the edges of the FOV, the limit is in- 
creased as expected. The faintest sources in the region have the 
PN-equivalent count rate about 1.5 x lO^'^ counts s^^, which we 
accept as our detection limit. In terms of flux, using the conver- 
sion above, this limit amounts to 3.3 x 10~^^ erg cm~^ s^^. 



3.3. Cross-identification 

We have cross-correlated the positions of the 235 sources with 
various optical and infrared catalogues. The 2MASS catalogue 
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Fig. 2. Bottom part of the distribution of the PN-equivalent count 
rates of the X-ray sources as a function of their distance d to 
the center of the FOV. The horizontal line shows the adopted 
detection limit. 



provides the most complete coverage of the Carina complex. We 
therefore selected this catalogue in order to determine the op- 
timal radius of cross-correlation. To this aim, we adopted the 
approach outlined by Jeffries et al. ( 119971 1 and applied to the 
XMM-Newton data of the young open clusters NGC 6530 and 
NGC 6383 by Rauw et al. (.20021 l2^3T l. The distribution of the 
cumulative number of catalogued sources as a function of the 
cross-correlation radius r is given by 



$(d <r) = A 



l-exp(^ 



{N- A) [l-cxp(-7rSr 



where N, A, a and B stand respectively for the total number of 
cross-correlated X-ray sources (A^ = 235), the number of true 
correlations, the uncertainty on the X-ray source position and 
the surface density of catalogue sources. For further details on 
the method we refer to the work of Jeffries et al. ( 1997). The fit- 
ting parameters A, a and B were obtained from the best fit to the 
actual distribution displayed in Fig.|3] For the 2MASS catalogue, 
we derive A ^ 217.6, a = 1.6"and B = 0.004 arcsec^^ j^e 
optimal radius that includes the majority of the true correlations 
while simultaneously limiting contamination by spurious corre- 
lations, is found to be around 4.0". We thus consider an infrared 
source as a possible counterpart if it falls within 4.0"of the coor- 
dinates of the X-ray source. This is significantly smaller than the 
corresponding optimal radius for the XMM data on NGC 6383. 
The main reason is the larger surface density of 2MASS sources 
in the Carina region. For r — 4.0", we expect statistically to 
achieve 208 true and 6 spurious correlations. 

We have also cross-correlated the positions of our X-ray 
sources with the AC03 XMM-Newton and EV03 Chandra X- 
ray surveys of the field, using a similar technique. The optimal 
cross-correlation radii were found to be around 5. 0"for the AC03 
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Fig. 3. Cumulative numbers of correlations between the X-ray 
detections and the 2MASS catalogue objects as a function of 
the correlation radius. The dotted curves correspond to the best 
fitting expressions for the real and spurious correlations. The 
dashed curve yields the sum of these terms and the dash-dotted 
vertical line corresponds to the adopted correlation radius of 
4.0". 



and 4.5"for the EV03 data statistically providing around 58/3 
and 54/5 true/false correlations respectively. 

The results of the correlations with the 2MASS and X- 
ray catalogues and the results of various optical studies of the 
Carina Nebula are listed in Table[3] The source designation in 
this table follows the naming conventions recommended by the 
lAU for serendipitous sources detected with XMM-Newton: the 
XMMU J prefix is followed by the right ascension HHMMSS.s 
(in hours, minutes, seconds and tenths of seconds, equinox 
J2000) and the dechnation of the source H-/-DDMMSS (in de- 
grees, arcminutes and arcseconds, equinox J2000), both trun- 
cated not rounded. The designations of the sources identified 
with the AC03 catalogue are taken from that catalogue; note that 
we still provide our own coordinates of such sources. 167 EPIC 
sources have a single 2MASS optical counterpart within a ra- 
dius of less than 4.0". The average angular separation between 
the X-ray source and the optical counterpart is (2.3 ± 0.9)". 
Cudworth et al. ( |1993l l performed a proper motion study of the 
open clusters in the Carina complex. Using photographic plates 
spanning about a century, they derived membership probabili- 
ties for 577 stars. We have cross-correlated the astrometry of our 
EPIC sources with the catalogue of Cudworth et al. Information 
on proper motion membership probabilities is provided in the 
last column of Table|3] 

Finally, given the limiting sensitivity of 1.5 x 
10^'^ counts s^^ for PN and the size of the field of view 
(roughly 0.25 square degree), we can use the logiV - log 5 
relation of Motch et al. (2003) for & ~ 0° to estimate the number 
of field stellar sources. In this way, we find that about 50 - 
70 field stars could be detected. Since the Motch et al. (2003]l 
relation was established from fields that do not harbour star 
formation regions, the significantly larger number of sources 



Table 3. Sample of the cross identification table of X-ray sources detected in our EPIC images of the Carina Nebula. 



x# 



XMMUJ 



RA 



[3] 



DEC 



[4] 



AC()3 EV03 



[5] 



[6] 



Nr. 
[7] 



[8] 



2MASS 
J H 



[10] 



Ks 

[11] 



V 

[12] 



B-V 

[13] 



U-B 

[14] 



Optical 
Name 

[15] 



Spectral Ref. Member 

Type (%) 

[16] [17] [18] 



1 


104236.8-594357 


10 42 36.83 


-59 43 57.9 


2 


104256.6-594220 


10 42 56.61 


-59 42 20.2 


3 


104303.2-593949 


10 43 03.21 


-59 39 49.3 


4 


104304.4-594902 


10 43 04.43 


-59 49 02.7 


5 


104307.8-593557 


10 43 07.82 


-59 35 57.4 


6 


104310.5-595012 


10 43 10.56 


-59 50 12.6 


7 


104310.9-595124 


10 43 10.91 


-59 5124.7 


8 


104311.4-594423 


10 43 11.40 


-59 44 24.0 


9 


104315.1-593344 


1043 15.12 


-59 33 44.2 


10 


104317.5-593321 


10 43 17.55 


-59 33 21.7 


11 


104319.8-594248 


10 43 19.86 


-59 42 48.5 


12 


104320.6-594619 


10 43 20.63 


-59 46 19.8 


13 


104325.9-594720 


10 43 25.93 


-59 47 20.5 


14 


104326.8-595251 


10 43 26.90 


-59 52 51.5 


15 


104327.4-594202 


10 43 27.48 


-59 42 02.4 


16 


104329.8-593344 


10 43 30.03 


-59 33 48.0 1 


17 


104333.7-594030 


10 43 33.76 


-59 40 30.0 


18 


104339.9-593445 


10 43 39.92 


-59 34 50.9 


19 


104340.7-594638 


10 43 40.75 


-59 46 38.4 


20 


104341.2-593355 


10 43 41.23 


-59 33 55.7 


21 


104341.3-593205 


10 43 41.49 


-59 32 09.8 3 


22 


104341.4-594223 


1043 41.61 


-59 42 26.8 4 


23 


104343.0-594411 


10 43 43.04 


-59 44 11.3 


24 


104343.1-593555 


10 43 43.16 


-59 35 55.7 


25 


104344.6-593459 


10 43 44.64 


-59 34 59.1 


26 


104345.0-595327 


10 43 45.08 


-59 53 27.7 


27 


104345.4-594159 


10 43 45.43 


-59 41 59.8 


28 


104345.6-593431 


10 43 45.60 


-59 34 31.6 


29 


104346.0-593221 


10 43 46.05 


-59 32 21.6 


30 


104346.3-593256 


10 43 46.36 


-59 32 56.6 


31 


104346.3-593116 


10 43 46.39 


-59 31 16.5 


32 


104347.3-593244 


10 43 47.35 


-59 32 44.2 


33 


104348.8-593324 


10 43 48.88 


-59 33 24.6 


34 


104349.2-593404 


10 43 49.22 


-59 34 05.0 


35 


104349.8-594453 


10 43 49.45 


-59 44 56.9 8 


36 


104351.0-595239 


10 43 51.03 


-59 52 40.0 


37 


104351.4-593117 


10 43 51.47 


-59 31 17.1 


38 


104351.6-593244 


10 43 51.63 


-59 32 44.3 


39 


104352.1-593556 


10 43 52.02 


-59 35 58.0 9 


40 


104352.3-593924 


10 43 52.48 


-59 39 24.7 10 


41 


104355.0-593130 


10 43 55.09 


-59 31 30.0 


42 


104355.2-593200 


10 43 55.23 


-59 32 00.7 


43 


104355.9-594926 


10 43 55.92 


-59 49 26.1 


44 


104356.0-593515 


10 43 56.03 


-59 35 15.1 


45 


104356.7-593252 


10 43 57.39 


-59 32 54.8 lie 
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3.9 



3.9 

3.7 

1.4 
3.3 
3.1 

2.9 
2.5 



3.0 

2.3 
2.3 

3.3 
3.5 

3.4 
3.9 

2.3 
3.3 
2.7 
0.1 
3.2 
2.6 
3.0 
1.4 
2.0 
2.0 
2.9 
2.4 
3.1 
1.9 
2.8 
1.7 
2.4 



13.00 12.37 12.18 



16.15 
13.76 

13.75 
8.70 

12.77 

15.15 
14.45 



13.41 
15.49 
13.21 

11.23 
12.70 

15.63 
14.41 

10.02 
14.37 
13.24 
15.26 

8.66 
12.37 
12.05 

9.55 
14.30 
12.67 
12.00 
13.56 
13.00 
13.28 
13.03 
13.26 
10.51 



14.41 
13.05 

12.97 

8.62 

12.30 

14.15 
13.52 



12.52 
14.69 
12.26 

10.80 
11.97 

14.31 
13.49 

9.99 
12.99 
12.42 
14.12 

8.52 
11.70 
11.56 

9.39 
13.13 
12.20 
11.60 
12.89 
12.27 
12.55 
11.91 
12.51 
10.04 



14.20 
12.83 

12.57 

8.50 

12.14 

13.58 
13.21 



12.21 
14.49 
11.96 

10.55 
11.71 

13.27 
13.12 

9.97 
12.18 
12.17 
13.51 

8.50 
11.51 
11.30 

9.26 
12.54 
12.06 
11.50 
12.63 
11.94 
12.22 
11.20 
12.26 

9.70 



12.47 



9.64 



15.43 



17.82 
15.72 



10.43 



9.65 

15.20 
10.73 

14.34 
13.75 

16.38 
16.58 

16.53 
12.52 



0.38 



0.26 



1.35 

0.75 
1.30 



0.14 



0.20 

1.36 
0.35 

0.77 
0.41 

1.26 
1.44 

1.46 
0.78 



-0.80 MJ3 



-0.76 HDE303316 



0.69 



MJ 



MJ 



DETWC-14 18 

DETWC-14 29 
DETWC-14 43 



DETWC 

DETWC 
DETWC 



DETWC-14 76, 65 
-0.77 MJ 126 

DETWC-14 78, 81 

-0.63 CPD-58°2611=Trl4 20 06V((f)) 



0.77 DETWC-14 107 
-0.55 Trl4 21 

Tr 14 Y 334 
-0.09 MJ 156 

DETWC-14 160 
DETWC-14 167 

DETWC-14 200 
0.43 Tr 14 14 



HD93129A + B 



09V 



The first and second columns yield the number of the X-ray source as well as the name according to the naming conventions for serendipitous XMM-Newton sources. To avoid duplicating the 
nomenclature, AC03 names are given for the sources cross-identified with the AC03 catalogue. The third and fourth columns provide the XMM-Netwon coordinates of the sources. Columns [5] and 
[6] list the star number in the X-ray catalogues of AC03 and EV03 according to our identification. Columns [7] to [11] summarize the results of the cross-correlation with the 2MASS catalogue. The 
columns labelled 'Nr' and 'd' yield the number of 2MASS counterparts within a 4.0" radius and the angular separation between the X-ray source and the IR counterpart. Columns [12] to [17] provide 
information on the properties of optical counterparts (when available). DETWC-14 and DETWC- 16 names yield the sequence number of the optical star in the catalogue of DeGioia-Eastwood et 
al. ( 1200 It for the Tr 14 and Tr 16 clusters respectively. In a similar way, Tr 14 and Tr 16 numbers refer to the numbering scheme introduced by Feinstein et al. ( |1973l l whilst Tr 14 Y and Tr 16 Y 
correspond to the convention of Cudworth et al. (1993V Finally, MI numbers are taken from Massey & lohnson tT993j . The column labelled 'Ref yields the reference for the optical photometry 
whilst the last column indicates the membership probability from the proper motion study of Cudworth et al. ( 119931 ). 
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detected in our Carina field clearly results from X-ray sources 
associated with the clusters in the field of view. 



3.4. Extragalactic background sources 

The line of sight towards i] Car (^n = 287.60°, bu = -0.63°) is 
nearly tangent to the Carina spiral arm. As a result, the neutral 
hydrogen column density along this direction must be quite large 
and should produce a substantial absorption of X-ray photons 
from extragalactic background sources (EBS). 

To get a rough first order estimate of the total Galactic ex- 
tinction along this line of sight, we made use of the DIRBE/IRAS 
extinction maps provided by Schlegel et al. ( I1998I I. Schlegel et 
al. caution that one has to be careful when using these maps 
near the Galactic plane. Moreover, the interstellar extinction in 
the Carina region is rather patchy and its properties could de- 
viate from those of the average Galactic extinction curve (e.g. 
Carraro et al. l2004l l. With these limitations in mind, we find that 
the DIRBE/IRAS maps indicate a rather large, but highly variable 
E{B — V) colour excess (between ^ 6.4 on average at the edges 
of the EPIC field and ~ 14.5 near 77 Car). Converting these val- 
ues into neutral hydrogen column densities by means of the gas 
to dust ratio of Bohlin et al. ( 119781 ). we estimate iVn in the range 
3.7-8.5 X 10^2 cm-2. 

Assuming that BBSs (most of which are probably Active 
Galactic Nuclei (AGN)) have a power-law spectrum with a pho- 
ton index of 1.4 (Giacconi et al. 12001). and suffer a total in- 
terstellar absorption of 3.7 - 8.5 x 10^^ cm^^, the above de- 
tection limits translate into unabsorbed fluxes of 1.0 - 1.5 x 
10"^'*ergcm"2s-i and 2.8 - 4.3 x 10~^'*ergcm"2 ^-i jjj jjje 

0.5 - 2.0keV and 2.0 - lOkeV band respectively. From the 
log A^ - log 5 relation of Giacconi et al. (2001), we expect to 
detect between ~ 16 and 24 extragalactic background sources. 



To discriminate these sources from the Galactic ones, the fol- 
lowing basic properties of AGNs could in principle be used (see 
e.g. Fiore . 19971 Nandra .2001. Bauer et al. .2004j : 



1. The spectra of AGNs are power laws with the canonical in- 
trinsic X-ray spectral slope F ^ 1.9 — 2.0. The average ob- 
served spectral slope of EBS is about 1.4 (Giacconi et al. 
I200TT ). 

2. Their flux is variable on timescales from a few hours to years. 



In practice, it is impossible to fit the individual source spectra 
since the candidates are very faint and the count statistics poor. 
For this reason, we have simulated the values of two PN hardness 
ratios HRi = (M - S)/{M + S) and HR2 = {H - M)/{H + 
M) for a grid of simple models. Within xspec, we have used 
the PN response matrices to generate fake spectra corresponding 
to differently absorbed power law spectra. We considered three 
photon indices F = 1.4, 1.7, 2.0 and several neutral hydrogen 
densities A^h = (0.1, 0.3, 1.0, 3.0, 7.0, 11.0, 15.0) x 10^2 cm'^. 
We then compared these simulated hardness ratios with the ob- 
served ones of all discrete sources not identified with optical or 
infra-red catalogues. For the list of the potential EBS candidates, 
we select those sources whose errorboxes cover the simulated 
hardness ratios. The resulting EBS candidates are the sources 3, 
6, 1 19, 124, 170, 185, 207. Their total number is seven which is 
lower than the estimation above; still the estimation gives only 
a rough idea about the number of the background sources, given 
the uncertainties involved. 



4. Open clusters in the Carina Nebuia 

The region covered by our EPIC data harbours several young 
open clusters that are extremely rich in hot massive stars; the 
most important ones being Trumpler 14 and 16. Tr 14 is a com- 
pact cluster containing three very hot 03 stars. Tr 14 is probably 
emerging from the eastern and near side of the dense parental 
molecular cloud (Tapia et al. 120031 ). Tr 14 might be slightly 
younger than Tr 16 as suggested by a fainter absolute magni- 
tude and a larger He II A 4686/ He II A 4541 ratio in the spectra 
of early and mid 0-type dwarfs (Walborn |1995l ). Collinder232 
is probably not a physical entity but a condensation of stars in 
the vicinity of Trl4 and 16 (see e.g. Walborn 1995 and Tapia 
et al. 2003). In a similar way, Collinder228 could actually be 
part of Tr 16. The apparent distinction between the two clusters 
might result from obscuration by a dust lane (Walborn 119951 ). 
Note an alternative view of Carraro et al. ( 120041 ) who argue 
that Collinder232 may be a physical aggregate (a rather sparse 
young open cluster), based on the comparison of its colour- 
magnitude diagrams and theoretical zero age main sequence 
(ZAMS) tracks. 

The extinction in the Carina complex has been subject to 
much controversy. Conflicting conclusions ranging from a nor- 
mal uniformly high to an anomalous and variable extinction have 
been proposed. This dilemma has important consequences for a 
consistent determination of the distance moduli (see e.g. the dis- 
cussion in Walborn 1995 ). In addition, since the line of sight is 
along the direction of a galactic spiral arm, there are probably 
many foreground and background objects that appear projected 
onto the Carina Nebula. 

The exact ages of the clusters and the star formation his- 
tory of the Carina Nebula are other tricky issues. Up until re- 
cent years, evidence for ongoing star formation in the Carina 
Nebula was rather scarce. Megeath et al. (11996 ) and Rathborne 
et al. ( 120021 ) reported evidence of star formation triggered in the 
molecular cloud exposed to the ionizing radiation of the hot stars 
of the Carina Nebula. None of the three near-IR sources N 1, N 3 
and N 4 proposed to be associated with embedded recently born 
O-stars is detected in our X-ray data. More recently. Smith et 
al. (2003) reported the discovery of dozens of candidates of so- 
called proplyds in the Carina Nebula, clearly demonstrating that 
the formation of low and intermediate mass stars is actively pro- 
ceeding. We note that none of the objects discussed by Smith et 
al. (2003 ) is detected as an X-ray source in our EPIC data. 

DeGioia-Eastwood et al. ( 120011 ) obtained UBV photometry 
of the Carina clusters. They found a significant population of 
pre-main sequence stars in Tr 14 and Tr 16 that have been form- 
ing over the last lOMyr whilst the most massive stars in Tr 16 
have ages between 1 and 3 Myr (with one exception, a star in 
the 5-6 Myr range in each cluster). According to their results, in- 
termediate mass star formation seem to have proceeded contin- 
uously over the past 10 Myr and the formation of intermediate 
mass stars started well before that of OB stars and was appar- 
ently not disrupted by the formation of these OB stars. Finally, 
Tapia et al. (120031 ) performed an extended UBVRIJHK pho- 
tometric study of the clusters in the Carina Nebula. They infer 
ages between less than 1 Myr and 6 Myr for Trl4 and 16. A 
small number of IR excess stars were found in both clusters. 
Tapia et al. further detected a population of 19 near-IR sources 
in the Car I Hll region whose formation could have been trig- 
gered by the action of the early-type stars. We note that none of 
these objects is detected in our EPIC data. 

Corcoran ( 119991 ) analysed ROSAT-HRI observations of the 
Carina complex. While he found the ix/Lboi ratio of early-type 
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stars to spread between 10^^ and 10^^, he did not find evidence 
for a population of X-ray bright pre-main sequence (PMS) stars. 
However, he cautioned that the failure to detect these objects 
could simply be due to the limited sensitivity of the ROSAT 
satellite. 

Trumpler 16 is particularly rich in massive binaries (see e.g. 
Levato et al. 1991) . The properties of several eclipsing or el- 
lipsoidal early-ty pe bin aries (HD 93205, Morrell e t al. ISOOTl 
Antokhina et al. |2000l Tr 16 -1, Freyhammer et al. |200T] and 
Tr 16-104, Rauw et al. I2001I I have been investigated recently. 
The components of the three binaries are all found to lie pretty 
close to the zero-age main sequence suggesting that they may 
be younger than 1 Myr (Rauw 120041) . The absolute magnitudes 
determined from the study of these eclipsing binaries suggest 
a distance modulus (DM) of 11.95 ± 0.06 (corresponding to a 
distance of 2.45 kpc) lower than the value of DM = 12.55 ± 
0.1 (which corresponds to a distance of 3.24 kpc) derived by 
Massey & Johnson (|1993l l and Massey, DeGioia-Eastwood, & 
Waterhouse (2001) from a spectoscopic and photometric study 
of the Trumpler 14 and 16 clusters. Allen & Hillier (1993]), 
Meaburn d 19991 1. and Davidson et al. (2001 ) investigated the ge- 
ometry and kinematics of the Homunculus Nebula and found 
the distance to 77 Car around 2.3± ~0.3kpc. Tapia et al. ( I2003I I 
derived the average distance to Trl4, Trl5 and Trl6 to be 
d — 2. 7 kpc from their UB Vi?l/HK photometric study, although 
they note a very large scatter in d. Carraro et al. ( 120041 ) confirm 
the distance of 2.5 kpc to Tr 14 but reach a different conclusion 
about the distance to Tr 16. Comparing the colour-magnitude 
diagrams and the Hertzsprung-Russel diagram (based on their 
[/BVi?I photometry) with post and pre-main sequence tracks and 
isochrones, they find that Tr 16 lies at the distance of about 4 kpc 
from the Sun. EV03 adopted the distance of 2. 5 kpc to the clus- 
ters in the Carina Nebula. In the current study we also adopt 
the distance of 2.5 kpc to all stars in the region. Note that while 
the absolute value of the X-ray and bolometric luminosities will 
change with the distance, their ratio remains unchanged. 

Fig. ID illustrates the JHKg colour-colour diagram of the 
2MASS infrared counterparts identified in our current study. We 
have included only those objects that were detected in all three 
filters J(1.25/im), iJ(1.65/xm) and Ks{2.n^m) with a pho- 
tometric accuracy better than a < 0.10 and which are single 
counterparts of EPIC sources within a correlation radius of 4.0". 
The intrinsic colours of dwarfs and giants (taken from Bessell & 
Brett |I98"8] l are indicated by the dashed and solid lines respec- 
tively. The open pentagons correspond to stars with J magni- 
tudes brighter than 10.0 (mostly moderately reddened early-type 
stars), whilst the crosses indicate fainter stars. 

In principle, this diagram can be used to identify pre-main 
sequence stars (Lada & Adams |I992| l either through a large 
circumstellar extinction characteristic of deeply embedded pro- 
tostars or through infrared excesses produced by circumstellar 
disks in classical TTauri stars (cTTs). The latter objects are ex- 
pected to fall mostly to the right of the reddening band. Meyer 
et al. (II997I I have shown that the locus of dereddened JHK 
colours of cTTs can be described by a simple linear relation be- 
tween J ~ H and H — K. This relation is shown in Fig.|4]by 
the dotted straight line. Only a modest number of the IR coun- 
terparts in Fig.|4]have colours that are consistent with these ob- 
jects being surrounded by large amounts of circumstellar mate- 
rial. The majority of the secondary X-ray sources could thus be 
weak line TTauri stars (wTTs) which display little or no evi- 
dence for circumstellar material in their spectra and photome- 
try. This situation is strongly reminiscent of other young clusters 
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Fig. 4. JHKs colour-colour diagram of the single IR counter- 
parts falling within the 4.0"correlation radius from an X-ray 
source. The heavy dashed and solid lines yield the intrinsic 
colours of main sequence and giant stars respectively, whereas 
the reddening vector is illustrated for Ay = 5 magnitudes of vi- 
sual extinction. Bright counterparts with J < 10.0 are displayed 
by open pentagons. The crosses indicate fainter IR sources. The 
dotted straight line yields the locus of dereddened colours of 
classical TTauri stars according to Meyer et al. ( |I9971 l. 



such as NGC6383 (Rauw et al. 12003] ) orNGC623I (Sanaet al. 
120071) where the X-ray selected pre-main sequence objects are 
also dominated by wTTs. 

5. Variability 

5.1. Short term variability 

To study possible variations of the detected sources within indi- 
vidual data sets, we applied the following variability tests: (i) we 
computed x^ to test the null hypothesis of a constant count rate 
level in the data; (ii) we applied the Kolmogorov-Smirnov (KS) 
test to check if the statistics of the count rate follows a normal 
distribution; (iii) we used the s,o-ca\\td probability of variability 
pov test suggested by Preibisch & Zinnecker ( 120021 ) in its modi- 
fied version suggested by Sana et al. C2004| l. 

In all methods, the source data were extracted from circu- 
lar regions with the radius equal to half of the distance to the 
nearest other source, or 60"if no nearby sources were present. 
The first two methods were applied to the background-corrected 
light curves of every source produced by binning photons in time 
and energy (see below). Note that, for faint sources, the KS test 
can sometimes fail for the background-subtracted light curves. 
The third method was applied on the series formed by photon 
arrival times (the photons were still binned in the same energy 
bins as for the x^ test). Thus, when using the pov method, it was 
not possible to account for the background. All variability tests 
were also applied to the background data. 

To use the x^ and KS criteria, we produced light curves of 
the sources with various time bin sizes (from 10 to 5000 sec) 
and in the soft S, medium M, hard H and total (0.4-lO.OkeV) 



10 



I.I. Antokhin et al.: XMM study of early type stars 



„ 0.100 

a 

»= 0.050 

o 0.000 
u 

S 0.100 

K 

- 0.050 

S 0.000 

0.060 
ffl 0.040 
0.020 
0.000 



K 



U 



o 0.015 
^ 0.010 

3 0.005 
o 0.000 

u 0.200 

^ 0.100 
3 0.050 

S 0.000 

u 0.200 
^ 0.150 
^ 0.100 

p 0.050 
S 0.000 

o 0.040 
„ 0.020 
g 0.000 

■^ -0.020 

ii 0.200 

CO 

□:; 

^ 0.100 

5 0.000 



r*^** ,\ \\ 



P+ti|fMlt^^^^'tMlli^' 



^Al ''^'^ 



|t.VM*/+f***+VV* ^*v*^ 



46 (3) 



,,^,,MiA+++^++**HtVW+++tHtty 



61 (4) 



f J| ?4 ( 1 ) 



# 



j*a\V,*V **^**^*WA^V^^AV^ 



80 (5) - 



'^*'**M.s** ,,,,,,,^ 



95 (5) 



88 (4) 



Ml 



J ^ , , , I , , , , I , , , , I ^ 

10" SxlO" 3x10* 

Time (sec) 



„ 0,300 

S 0.200 
c 0.100 
S 0.000 
u 0.060 
ffi 0.040 
c 0.020 

3 0.000 

S 0.080 
" 0.060 
I 0.040 
g 0.020 

a; 

o 0.040 

^ 0.020 

c 

S 0.000 

„ 0.030 

E 0.020 

§ 0.010 

" 0.000 

i! 0.060 

^ 0.040 

§ 0.020 

S 0.000 

aj 

o 0.040 

.^ 0.020 

o 0.000 

u 

■g 0.100 
J 0.050 
D 0.000 



'"'■' K, 



\*/Jx^*\^ii^*,f^4^M\**} *W 






145 (1) 



••."• **s 






189 (5) I 



:^+^n^^+ ++t+ttt,+Mt^^^f++^^^l 






+++tH^'t ^ 



194 (1) ^, i + i| 

^ — ^ — \ — \ — \ — \ — \ — \ — \ — \ — \ — \ — \ — \ — \ — \ — 



195 (4) 



10" 2x10" 3x10" 

Time (sec) 



Fig. 5. PN background-subtracted light curves in counts s^^ of sources showing significant short term variability (MOS 1 light curve 
for the source 73). The energy range is 0.4-lO.OkeV. The source numbers followed by the data set numbers in parantheses are 
shown (for source identification see text and Table|3]l- In this and the following light curve plots (Fig|6]-[8]l the time shown is relative 
to the start of the exposure. 



energy ranges. Since the Carina region is quite crowded, it was 
often impossible to find a source-free area in the vicinity of a 
given source to extract a background light curve. For this rea- 
son, we adopted the following scheme when extracting the back- 
ground light curves. For the MOS instruments, the internal in- 
strumental background is more or less constant within a single 
CCD. For the PN, the internal background primarily varies with 
the distance along the detector "y"-axis from the central detec- 
tor "x"-axis (separating 2 sets of the PN CCDs). Thus, for the 
MOS instruments, we defined several object-free circular areas 
on every CCD and extracted one background light curve per en- 
ergy range per MOS CCD from all these areas. The background 
light curves from a given MOS CCD were used for all objects 
located on that CCD. We performed a similar procedure with 
the PN background light curves except that circular background 
areas were selected in 10 strips parallel to the central PN "detx"- 
axis. Apart from instrumental background, there exists diffuse 
X-ray emission from the Carina nebula. Thus the background 
light curves obtained as above, may overestimate or underes- 
timate the real background level near particular sources. This 
could influence the results for the faintest sources. However, in 
practice this is not very important as (i)we use several meth- 
ods for detecting variability, including those without background 
subtraction; (ii) the count rates of the faintest sources are small 
anyway resulting in large errors and small chances to find vari- 
ability; (iii) as the diffuse emission has large spatial scale (about 



the size of one CCD or larger) the adopted procedure roughly 
accounts for global changes of the diffuse background. 

To increase signal, the sum of the background-corrected 
MOSl and MOS2 light curves was also analysed whenever pos- 
sible. In Fig|5]we show the PN light curves of those objects for 
which the null hypothesis is rejected at a confidence level of 99% 
using all three variability test methods (for source 8 1 we show 
the MOS 1 light curve as the source falls onto a gap between PN 
CCDs). We do not show the MOSl and MOS2 light curves due 
to lack of space. The behaviour of these light curves is similar 
to the PN ones. The time bin sizes of the plots were selected ac- 
cording to typical variability time scales of the corresponding 
light curves. 

Little is known about the nature of most of the variable 
sources. Only one source (#189) is identified with HD 93343 
(07V(n)). The sources #36 (=MJ 156), #46 (=Trl418), #104 
(=DETWC-16 10,5), #145 (=Tr 14 Y66), and #195 (=DETWC- 
16 345) are identified with the 2MASS and optical catalogues 
but no spectral type is known. The sources #74, 80, 115, 137, 
177, 194 are identified with the 2MASS catalogue, no optical 
identification is known. The sources #61, 73, 88, 95, 193 are not 
identified with any optical or IR catalogues. All the sources dis- 
played show variability at different levels, and 3 of them (#74, 
80 and 104) are clearly flaring. 
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Fig. 6. PN light curves of source #74 in the soft, medium, and Fig. 8. PN light curves of source #104 = DETWC-16 10,5 (data 
hard energy bands. set 1) in the soft, medium, and hard energy bands. 
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Fig. 7. PN light curves of source #80 in the soft, medium, and 
hard energy bands. 



5.1 .1 . The probable pre-main sequence stars #74, #80, 
and #104. 

Further insight into the nature of the flare-type sources may be 
obtained from the analysis of their light curves and spectra dur- 
ing the flare events. But first, let us consider the IR properties of 
the flaring X-ray sources. 

Fig.|9]illustrates the K \s. J — K colour-magnitude diagram 
of the 2MASS counterparts in the EPIC field of view. We used 
the March 2003 update of the colour transformations, initially 
derived by Carpenter ( 1200 11 1 and available on the 2MASS web- 
sitqj, to convert the 2MASS colours and magnitudes into the 
homogenized JHK photometric system introduced by Bessell 
& Brett (1988). 

Assuming all the sources are at the same distance, the 
sources to the right of the main sequence could indicate the ex- 
istence of a population of pre-main sequence stars. The counter- 
parts of three flaring X-ray sources are illustrated by the triangles 
(downward triangles stand for the two possible counterparts of 
source #104, whilst the upward triangles show sources #74 and 
#80). The flaring X-ray sources with single 2MASS counterparts 
are clearly located at positions consistent with these objects be- 
ing PMS stars. 

The X-ray light curves of the three flaring sources are shown 
in Fig|6]-[8] It is immediately evident that the flares occur mostly 
in the hard energy range which is consistent with these observed 
flares being produced by the coronal activity in PMS stars. 

Assuming that the plasma producing the observed variabiltiy 
is confined in a closed coronal loop as is the case for solar flares, 
we can use the method proposed by Serio et al. ( 119910 . These 
authors established an analytical relation between the loop half- 
length I, the maximum temperature at the top of the loop Tmax 
and the thermodynamic decay time r (see also Briggs & Pye 
120031 Giardino et al. |2004l and Favata |2505T l. 

We have determined the 1/e decay times of the flares of 
sources #74, 80 and 104 by a least square fit of the exponen- 



|http : //www. ipac ■ caltech . edu/2mass/ index . html| 
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Table 4. Properties of the flaring sources. 



Source 


Rev. 


T 


1^-'- obs 


I X F(C) 


fx (0.4-lOkeV) 


J - K 


K 






(ksec) 


(keV) 


(R©) 


(ergcm-^s-i) 






74 


115 


1.5 ±0.3 


13.3 ±7.5 


1.1 ±0.3 


9.3 X 10-" 


1.452 


14.198 


80 


285 


2.5 ±0.6 


4.2 ±1.1 


0.9 ±0.2 


2.9 X 10"" 


1.381 


12.718 


104 


115 


2.1 ±0.4 


4.6 ±0.7 


0.8 ±0.1 


9.1 X 10"" 







r yields the 1/e decay time, fcTobs is the best fit plasma temperature during the flare and I is the loop half-length. The X-ray fluxes correspond to 
the spectrum during the flare. 
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Fig. 9. K vs. J — K colour-magnitude diagram of the 2MASS 
counterparts with good quality near-IR photometry of the X-ray 
sources in the EPIC field of view. The open pentagons indi- 
cate early-type stars, whilst the triangles stand for flaring X-ray 
sources. The locus of the un-reddened main sequence is indi- 
cated for a distance modulus of 12.0 by the leftmost solid line. 
The reddened main sequence is indicated for two different red- 
dening values Ay — 5 (middle solid line) and Ay = 10 (right- 
most solid line). 



tial decay phase of the observed light curves. Considering the 
typical count rates of the sources studied here, the most use- 
ful time bins to study the light curves are 500 and 200 s for the 
MOS and PN data respectively. The decay times derived from 
the fits of the different EPIC light curves agree within 20%. 
The results are given in Table|4] All three sources have rather 
short decay times indicating that the X-ray emitting loops must 
be relatively compact. We fitted the spectra of the sources dur- 
ing the flare with an absorbed single temperature apec model. 
The best-fit column densities of the three sources are found to 
be 0.26 - 0.39 x 10^^ cm"^. The best-fit temperatures are quite 
high (kT > 4.2 keV) as expected for pre-main sequence stars. 
The measured decay time of the flare actually exceeds the in- 
trinsic thermodynamical decay time due to the effect of heating 
that can extend into the decay phase of the flare. The relation be- 
tween the two F{() — Tobs/Tth can in principle be determined 
from the slope ( of the flare decay in a log T - log \/EM di- 
agram (Favata et al. I2005I I. Unfortunately, the flaring sources 



discussed here were too faint to allow a time-resolved spectral 
analysis. Therefore, we can only determine I x F{C,). For flar- 
ing loops that are freely decaying with no heating, F{Q ~ 2 
(Favata et al. 2005) which actually represents the lower limit on 
F{C). Therefore, the typical loop sizes of the three sources are 
<O.55R0(seeTable|lli. 



5.2. Set-to-set variability 

To detect possible variability of the X-ray sources between the 
data sets we performed a x^ tsst at three different confidence 
levels (90, 95 and 99%) to test the null hypothesis of a constant 
count rate level in the data. The tests were performed for the 
count rates in the soft S, medium M, hard H, and total (0.4- 
lO.OkeV) energy bands as well as for the hardness ratios HRi 
and HR2, for all EPIC instruments. The results of the test are 
included in Table|2] (second column). The variability status is 
assigned as follows. If the null hypothesis is rejected at a con- 
fidence level of 95% in the total energy range in all instruments 
simultaneously, the source is labeled "Var". If the null hypoth- 
esis cannot be rejected in all instruments, the source is labeled 
"Const". If the null hypothesis is rejected in some instruments 
only, the variability status is set to "Uncert". Finally, if there is 
not sufficient data to perform the test, the variability status is set 
to "Nolnfo". 

Of course, given the limited number of data sets, the detected 
variabiltiy of a source tells us hardly more than the mere fact 
that the source is variable on a time scale of one day or one 
year (the time intervals between our data sets). The nature of 
this variability for the majority of the sources remains uncertain. 
The purpose of this test was to set a preliminary base for future 
studies. Some sources like HD 93205 (#1 17) are early-type bina- 
ries, in which case the X-ray flux varies with the orbital motion. 
Spectral analysis may help in disentangling the nature of the de- 
tected sources and their variability. Such an analysis will be the 
subject of forthcoming papers. 



6. X-ray properties of early-type stars 

We extracted the spectra of all early-type stars detected in the 
images, using the same apertures as for their light curves. The 
background spectra were extracted from circular apertures cho- 
sen in source-free areas as close as possible to the sources under 
investigation, so that the diffuse background emission is simi- 
lar to that at a source position. The corresponding redistribution 
matrices and ancillary response files were generated for every 
source, the spectra were rebinned so that they contain a mini- 
mum of 10 counts per channel. 
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x# 



Source Name 



Spectral Type 



E{B-V) A^H^i Nn kT 

(10^2 cm-2) (io22 cm-2) (keV) 



Norm 

(10-4) 



xVdof 



dof 



MJ126 


B2Vb 


0.42 


0.24 


CPD-58°2611 


06V((f)) 


0.45 


0.26 


Trl4 21 


09V 


0.61 


0.35 


HD93129A+B 


03Iab...+03.5V 


0.56 


0.32 



0.24 



46 


Tr 14 18t 


51 


HD 93130 


59 


Tr 16 124 


63 


HD 93160 



BOV 
06III(t) 

BIV 
06III(t) 



67 HD93161A+B (08 V+09 V)+06.5 V(f) 

71 HD93162(WR25)tt WN6h 



108 CPD-58°2649 



110 



117 



142 



178 



HD 93204 



HD 93205 



126 HDE303311 



135 CPD-59°2600 



HD 93250 



147 CPD-59°2603 



CPD-59°2626 



179 HDE303308ttt 



183 


CPD-59°2629 


189 


HD 93343 


191 


CPD-59°2636 


197 


CPD-59°2641 


218 


MJ596 


235 


HDE 303304 



BOV 

05V((f)) 

03.5V+08V 

05V 
06V((f)) 

03V 



03V((f)) 



08.5 V 

07V(n) 

07 V+08 V+09 V 

06V((f))+? 
05 V((f))+09.5 V 

0+... 



0.52 
0.48 
0.50 
0.31 

0.29 

0.69 

0.54 
0.39 
0.40 
0.44 
0.48 
0.48 



07 V+09.5 V+B0.2 IV 0.43 

0.43 

07V 0.65 



0.44 

0.78 

0.54 

0.56 
0.58 
0.73 

0.61 



0.30 
0.28 
0.29 
0.18 

0.17 

0.40 

0.31 
0.23 
0.23 
0.26 
0.29 
0.29 



0.25 
0.25 
0.38 

0.26 



0.45 

0.31 

0.32 
0.34 
0.42 

0.35 



-0.06 
-0.00 

"■'"-0.12 
1 91-1-0.14 

"■^"-0.00 
"■' '-0.06 



0.16^ 
0.68^ 



R+0.06 
-0.03 
fO.22 
-0.17 

^••^^-0.00 

0.28+«-°2 



0.00 

+0.06 

0.00 

+0.09 

0.07 

+0.39 

0.35 

0.08 

0.06 

61+0-57 
^•"r_Q 32 



0.29 
0.58 
0.35 
0.57 



+0.09 
0.09 

+0.13 
■0.10 

+0.06 
-0.00 



0.00^ 



0.59 
0.29 
0.00 

0.78t«:li 

fO.44 
-0.00 

34+0-18 
'J-^^-o.ii 

0.42 
0.03 
0.26 



+0.04 
0.15 
0.13 

-0.01 
0.09 
0.00 

00+0-2* 
'^•'-"J-0.00 
,+0.02 
-0.00 
,+2.51 
-0.00 
0.02 
0.00 
Q n, +0.06 
'^••^^-0.05 



0.29 
0.00 
0.29 



0.00 
0.26 



+0.03 
-0.00 

+0.21 
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0.47 
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1.01 
0.16 
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-0.04 

+0.08 
-0.26 

+0.22 
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6.82 

0.37 



-0.06 
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^••-'^-0.32 

-1.17 
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+0.36 

-0.46 

+0.82 

-0.32 
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"••-'^-0.20 

3+6.60 
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3+0.95 

-1.00 

-|p;+0.02 

^•l'^-0.02 

+3.43 

2.82 

+0.27 

0.23 

+ 1.61 

1.14 

+0.36 

0.08 

+0.39 

0.25 

+0.23 

0.23 

+0.69 

0.37 

+ 1.02 

0.40 

17.4 

10.1 



1.19^ 



171: 

4.79"+ 



9.18 
3.16 
2.45 
0.47 
1.37 
0.95 
0.54 
1.13 
11.0 



1.14 


135 


1.69 


51 


1.45 


225 



1.11 



1.05 



1.05 



1.02 



I.IS 



1.13 



1.05 



1.11 



1.01 



1191 



1.22 


185 


1.13 


263 


1.16 


101 


0.99 


258 



342 



1017 



92 



239 



524 



238 



275 



1408 



1.20 165 

1.11 55 
1.13 632 

0.96 924 

1.33 159 

1.30 120 

1.12 102 

1.21 160 

0.96 107 



^ISM '^ ''1^ ISM column density of neutral hydrogen. Norm is the normalization factor of a thermal model component given by Norm 

where d is the distance to the source (cm), Ue and nu are the electron and H densities (cm^ ^). 

t - Parameters may be unreliable as the star is located in an area of extremely varying background and close to a cluster of bright stars. 

tf -Metal abundance 0.59±0.07. 

t f f - Parameters may be unreliable due to proximity to rj Car and a strong non-uniform background. 
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6.1. Spectral models 

Table|5]lists all detected early-type stars along with the parame- 
ters of the spectral fits. Spectral fitting was generally performed 
for all data sets and all instruments MOSl, MOs2, and PN simul- 
taneously. In some cases when a source was variable or the data 
from some instruments in a particular data set were affected by 
the CCD gaps, a subset of the data was used. To quantify the 
physical properties of the X-ray emission we adjusted a series 
of optically thin thermal plasma apec models to the observed 
spectra using xspec. The general form of a multi-temperature 
model was wabsi* (apeci+wabs2*apec2+ . . . ) ■ wabsi is 
an absorption model using the cross-sections for cosmic abun- 
dances. We required that the column density of neutral hydrogen 
for this model component was not smaller than the value given 
in the 5th column of Table|5] The column density was calcu- 
lated individually for every star using the empirical relation from 
Bohlin et al. (1978). Additional wabs components account for 
possible internal absorption of X-ray radiation in the winds of 
early-type stars. Metal abundances were fixed to the solar value 
in all cases except WR25, in which case it was a free parame- 
ter. In a few cases an Fe line complex (or indication of it) was 
present in the spectra at 6.4-6.8 keV. This line complex may in- 
clude the Fe 1 fluorescence line at ^6.4 keV and the lines of Fe 
XXV and XXVI centered around 6.7 keV The potential pres- 
ence of the Fe 1 line (indicative of cool material) may skew the 
high temperature model component required to fit the Fe XXV 
and XXVI lines. We did a check for the presence of the Fe I line 
as follows. We start fitting with just one thermal model compo- 
nent and include a gaussian component at 6.4 keV If necessary, 
we increase the number of thermal components until a satisfac- 
tory fit is achieved. We then remove the gaussian component and 
check if the goodness of fit decreases significantly. If not, we per- 
manently remove this component from the model. It turned out 
that no significant Fe I emission was found in any sources. 

Table|5] does not include the brightest object in the FOV, 
rj Car. The nature of this star is so complicated that it makes little 
sense to give the results of a necessarily superficial analysis here. 
A detailed study of XMM-Newton observations of this star was 
recently published by Hamaguchi et al. (2007) . A few weak OB 
stars were also excluded from the table as no reliable spectral fits 
could be performed. 

The number of thermal components for different sources 
varies from one to three. Spectra of some stars required rather 
high temperature model components. Below we will discuss pos- 
sible reasons for this. We also tried to fit the spectra with com- 
binations of thermal and power law models. These did not result 
in better fits. Note that for some weak sources, the fitted param- 
eters are not really well constrained and about the same quality 
of the fit can be obtained in a single-temperature model for dif- 
ferent combinations of Nh and kT. However, as we are most 
interested in the source ^wxes within the XMM-Newton energy 
range, the somewhat formal fits are acceptable. Our model pa- 
rameters are roughly in agreement with those obtained by AC03 
and EV04 (the latter are based on a Chandra observation of the 
region). The exact comparison is impossible as we used some- 
what different models (e.g. different number of temperatures in 
thermal models), owing to different quality of the data. 

A few examples of the observed spectra along with the fitted 
models are shown in Fig [T0l - [T3] These examples show WR25 
and some OB stars of different brightness, to illustrate the quality 
of the data. 




channel energy (keV) 

Fig. 10. EPIC spectra of WR 25 from the data set 3 along with 
the best-fit model. In the top panel, the upper and lower data 
correspond to the PN and MOS spectra, respectively, and the solid 
lines yields the best-fitting model. The bottom panel shows the 
contribution of individual energy bins to the x^ of the fit. 




1 2 

channel energy (keV) 



Fig. 11. Example of a weak source. PN spectra of the source #108 
(CPD-58°2649) from the data sets 1,2, and 5. 




1 2 

channel energy (keV) 



Fig. 12. PN spectrum of the source #117 (HD 93205) from the 
data sets 4 and 5. 
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Fig. 14. Lx versus Lboi for early-type stars in the Carina Nebula in different energy bands. B-type stars and WR 25 are excluded 
from the least-square fits. The dashed lines show the best-fit approximations of the Lx - iboi relation for O (dashed lines) and B 
(dotted Unes) stars in NGC 6231 (Sana et al. l2006t (hereafter S06). 
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channel energy (keV) 



Fig. 13. PN specti-um of the source #142 (HD 93250) from the 
data sets 3-5. Note the presence of the Fe line complex at 6.4- 
6.7 keV. 



6.2. X-ray fluxes and L^/Lhoi ratio 



The X-ray absorbed and unabsorbed (i.e. corrected for the in- 
terstellar absorption) fluxes and bolometric luminosities of OB 
stars detected in the XMM-Newton data are presented in Table|6] 
The first column is the source number, the second to fifth 
columns are the absorbed X-ray fluxes in the soft, medium, hard, 
and total energy bands, while the following four columns are the 
unabsorbed fluxes in the same bands. The tenth column is the 
bolometric magnitude. 

EV03 thoroughly discuss bolometric magnitudes for the 
early type stars they detected in their Chandra data. We adopted 
their bolometric magnitudes for the sources in common. For the 
stars not detected in the Chandra data, we, like EV03, com- 
pute the bolometric magnitudes using the approach and formu- 
lae from Massey et al. ( |2000| l. Briefly, the reddening Ay = 
R X E{B — V) with R = 3.2 was computed for each star 
to get its de-reddened visual magnitude. The bolometric cor- 
rection of a star was computed using the relationship BC — 
27.66 - 6.84 x log Toff- In the complex case of HD93161 con- 
sisting of three components Aa, Ab, and B, we computed the 
bolometric luminosities of every component separately using the 
information about the colors and temperatures from Naze et al. 
(2005) and then summed them up to get the total bolometric 
luminosity of the system. For MJ 596 we used the bolometric 
magnitude obtained by Niemela et al. ( |20061 l from the detailed 
analysis of spectral and photometric data. 

Fig. [14] shows the intrinsic Lx versus Lboi for the stars 
in Table|6l in every energy band. It is known from the litera- 
ture (e.g. S06) that O- and B-type stars show a different be- 
haviour. For this reason we marked them with different sym- 
bols. Presumed single O-type stars are marked with solid cir- 
cles, binary stars are plotted as solid triangles, and B-type stars 
are marked with open circles. The demarcation line between two 
types of stars is located at about Lboi — 10'^* erg s^^. We shall 
not discuss B-type stars here as their number in our sample is 
small. 

The X-ray and bolometric luminosities of O-type stars are 
clearly correlated in the soft and medium bands. In the hard band 
the errors are large so no definite conclusion can be drawn. 

The solid lines in Fig.[T4lrepresent the least-square fits of the 
observed distributions of O-type stars (both binary and single) in 
the form 



log(Lx) = a + log(Lboi) 

which represents a simple scaling relation Lx = lO'' x Lboi- 
Our data do not allow us to search for a more sophisticated rela- 
tion e.g. in the form of a power law. 

As seen from Fig. [14] there is no clear distinction between 
binary and single O-type stars in the ix-^boi plane. We per- 
formed various tests on the two samples taken separately. It 
turned out that the correlation coefficients for the two sam- 
ples, the parameters of the least square fits, and the fit residuals 
were all very close. The numerical values of the least-square un- 
weighted fit parameters for the whole (single plus binary) sample 
of O-type stars are: 



log(L 



soft\ 



iogiLT") 






(-6.82 ± 0.56) + log(Lboi) 
(-7.18 ± 0.92) + log(Lboi) 
(-6.58 ± 0.79) + log(Lboi) 



(1) 



The errors of the parameter are obtained assuming a good fit 
(Press et al. 1 19921 1 and are rather formal. The dispersion of the 
luminosities around Eqs. 1 are 29%, 48%, and 41% in the three 
bands above. 

In the soft and medium bands, a few O-type stars devi- 
ate from the least-square scaling relation more than the others. 
These stars are labeled with their names in Fig. [14] There is no 
unique reason for this deviation. HDE 303304 shows increased 
X-ray luminosity in the soft band. The star is somewhat unusual; 
it is classified as "OBH-(le)" (that is having at least one emission 
line in its spectrum) by Stephenson & Sanduleak ( |1971| l. It is 
not clear if this can be related to the increased luminosity in the 
soft band. On the other hand, the binary hypothesis is not very 
likely. If the increased X-ray flux is due to wind-wind collision, 
then the X-ray spectrum should also be harder than a typical 
spectrum of a single star, and thus the X-ray limunosities in the 
medium and hard X-ray bands ought to be higher than average. 
Evidently this is not the case. More likely is that the soft X-ray 
flux of HDE 303304 is affected by the diffuse background X-ray 
emission. The source is weak and the diffuse X-ray emission has 
maximal flux in the soft X-ray band. It is also very inhomoge- 
neous across the field of view. 

The 08.5 V star CPD -59°2629 (=Tr 16 22) is also discov- 
ered as an unusually bright and hard X-ray source for its spec- 
tral type by Chandra instruments (EV04). The latter authors sug- 
gested that the unusually large Lx/L\^oh as well as a surprisingly 
hard spectrum, could indicate the presence of a wind confined 
by the magnetic fields of this star Another possible explanation 
could be the presence of a low-mass pre-main sequence compan- 
ion. These objects can have Lx/iboi up to 0.001 and even 0.01 
(during flares). However, in the absolute terms, it is unlikely that 
they would be bright enough to be seen against the intrinsic X- 
ray emission of an early-type star E.g., InNGC6231, PMS stars 
can typically display Lx up to 10^^ erg s~^ (Sana et al . 120071 1. 

The deviations in the soft and medium bands of Tr 1421 is 
likely related to the fact that this weak source is located near 
the edge of the field of view, in an extremely crowded area with 
strong and variable background diffuse emission. 

Despite the fact that some deviations from the least-square 
scaling relations in Fig. [T4| could be understood as instrumen- 
tal effects we did not exclude the corresponding stars from the 
fits. The reason is that, even for stars which are visually "well- 
behaved" , in the plots the uncertainies in the exact value of the 
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Table 6. X-ray fluxes and absolute magnitudes of early type stars in the Carina Nebula. 
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x# 




Absorbed Flux 






Unabsorbed Flux 




Mboi 




[0.4-1.0] 


[1.0-2.5] 


[2.5-10.0] 


[0.4-10.0] 


[0.4-1.0] 


[1.0-2.5] 


[2.5-10.0] 


[0.4-10.0] 




26 


3.52 


2.43 


1.37 


6.09 


9.12 


3.40 


0.142 


12.66 


-6.21 


30 


0.86 


1.74 


0.10 


2.71 


2.22 


2.41 


0.104 


4.74 


-6.30 


33 


1.28 


7.22 


8.65 


17.1 


4.88 


9.94 


8.92 


23.72 


-5.95 


45 


32.4 


74.5 


33.5 


140.4 


151.5 


105.8 


34.6 


291.8 


-11.5 


46 


1.80 


11.0 


11.0 


23.8 


8.57 


14.81 


11.30 


34.7 


-3.86 


51 


4.00 


2.66 


0.11 


6.77 


12.7 


3.99 


0.12 


16.8 


-8.83 


59 


0.369 


2.18 


1.98 


4.53 


1.68 


2.92 


2.03 


6.63 


-5.76 


63 


4.91 


4.58 


0.94 


10.4 


11.0 


5.75 


0.96 


17.7 


-9.4 


67 


6.71 


5.91 


0.97 


13.6 


14.0 


7.36 


0.99 


22.4 


-9.39 


71 


46.9 


216.1 


150.0 


413.0 


205.0 


324.9 


155.5 


685.4 


-10.1 


108 


1.23 


1.18 


1.07 


3.48 


6.40 


1.67 


1.10 


9.17 


-6.9 


110 


5.99 


3.29 


0.16 


9.43 


20.6 


4.51 


0.16 


25.2 


-9.0 


117 


18.2 


11.4 


0.52 


30.0 


58.1 


15.5 


0.55 


74.2 


-10.1 


126 


3.34 


2.26 


0.84 


6.44 


13.4 


3.17 


0.86 


17.4 


-8.6 


135 


5.90 


4.76 


1.66 


12.3 


19.4 


7.04 


1.69 


28.1 


-9.0 


142 


32.8 


59.0 


44.0 


135.8 


124.0 


81.5 


45.0 


250.5 


-10.7 


147 


2.02 


1.45 


0.10 


3.55 


6.60 


2.02 


0.10 


8.72 


-8.5 


178 


1.62 


2.35 


0.40 


4.36 


11.3 


3.60 


0.42 


15.3 


-7.5 


179 


15.9 


17.7 


89.6 


123.2 


70.6 


23.0 


90.9 


184.5 


-9.8 


183 


3.00 


12.4 


9.88 


25.3 


18.2 


19.8 


10.3 


48.2 


-6.8 


189 


1.97 


3.03 


2.57 


7.54 


7.48 


4.36 


2.63 


14.3 


-8.0 


191 


3.39 


3.87 


0.20 


7.46 


11.6 


5.99 


0.21 


17.8 


-8.3 


197 


2.23 


3.02 


0.22 


5.47 


7.71 


4.68 


0.23 


12.6 


-8.7 


218 


1.27 


3.63 


5.41 


10.3 


6.47 


5.45 


5.58 


17.5 


-8.18 


235 


4.53 


1.08 


0.001 


5.61 


35.2 


1.88 


0.0012 


37.1 


-7.45 



Fluxes are in units 10 erg cm 



diffuse background emission, and especially in the inferred neu- 
tral hydrogen column densities, may still be important factors 
influencing the exact location of the data points in Fig. [141 We 
believe that excluding the deviating points for the sole purpose 
of getting a nicer- looking plot is not a good approach in this case. 



6.3. Comparison with X-ray fluxes from ottier reviews 

Before we turn to the discussion of the results let us briefly check 
the consistency of the fluxes obtained in the current study, with 
the previous values presented in EV03 and AC03. 

First, note that five OB stars from the list of EV03 (#8, 9, 39, 
44, 141) are not detected in our XMM-Newton data, even de- 
spite the fact that the effective telescope area and exposure times 
are larger in the XMM case. This is explained by the lower spa- 
tial resolution of XMM. All these stars are weak X-ray sources. 
Stars #8 and #9 are located near the brightest object in the region 
- T] Car (about 39"and 57"respectively), and are lost in the bright 
background created by the Homunculus. The EV03 source #44 
is between the XMM close sources #189 and #191, the EV03 
source #141 is very weak and is located near another XMM 
source #90. As a consequence, these sources were not detected 
by the SAS detection routine. 



A comparison of absorbed XMM and Chandra fluxes for 
common early-type sources is shown in Fig. [15] EV03 com- 
puted their unabsorbed Chandra fluxes in the energy band 0.5- 
2.04 keV from the observed count rates assuming a single- 
temperature plasma model with the temperature kT=0.385keV 
and the hydrogen column density Nh — 0.3 x 10^^ cm~^. 
We converted these fluxes (taken from the Table 3 of EV03) to 
the absorbed ones using the above model parameters. XMM- 
Newton absorbed fluxes were computed in the same energy band 
0. 5-2.04 keV using the best-fit models from Table|5] The results 
of the two telescopes are in quite good agreement, especially 
accounting for the calibration issues at the time of Chandra ob- 
servation. 



A comparison of absorbed fluxes of early-type stars from the 
current study with the fluxes from AC03 (their Table 3) is shown 
in Fig. [16] Two weak stars with problematic fluxes (our #179 and 
#197) are excluded from the plot. The AC03 fluxes are computed 
in the energy range 0.3-12keV which explains the systematic 
shift between our and their fluxes. Otherwise, the agreement is 
quite good. 
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7. Discussion 

7. 1. General properties of the field 

The X-ray images of the Carina region are dominated by dif- 
fuse emission and strong discrete sources. The brightest discrete 
sources are associated with the hot massive stars in the various 
clusters that populate this region (see Sect.|4]i. Among the fainter 
sources, there are likely some extragalactic background objects, 
but the vast majority are likely to be low-mass stars, either fore- 
ground coronal sources or pre-main sequence stars in the Carina 
complex (see Sect. 15.1. Ti l. Compared to some other open clusters 
such as NGC 6231, there are less PMS candidates detected. This 
could be due to a lower detection efficiency (due to the strong 
diffuse emission) and the larger distance of the Carina complex 
compared to NGC 623 1 . 



-14 
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Fig. 15. Comparison of absorbed XMM-Newton and Chandra 
fluxes in the energy band 0.5-2. 04 keV. 



7.2. Early-type stars 
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Fig. 16. Comparison of absorbed fluxes from the current study 
(energy band 0.4 - lO.OkeV) and the study of AC03 (energy 
band 0.3 -12.0keV). 



In Section 16.11 we derived the spectral characteristics and the 
observed ix-^boi relations for 0-type stars in our sample. 
Recently, a similar investigation was pulished regarding the 
early-type population in NGC 6231 by S06. The authors also 
provide a review of the previous works on this subject. The 
S06 linear fits of the Lx - iboi relations for O and B stars in 
NGC 6231 are shown in Fig. [14] by dashed and dotted Unes re- 
spectively. 

When comparing our results with those of S06 one should 
keep in mind that the Carina Nebula is known as a region with 
very patchy interstellar absorption, thanks to numerous arms of 
dust and gas penetrating through the volume of the stellar asso- 
ciations. For this reason it is extremely difficult to account accu- 
rately for the ISM absorption. The situation with NGC 623 1 is 
less complex and this is a probable explanation for the scatter of 
data points in our Fig. [14] which is somewhat larger than in the 
case of SOqJ- In addition, the difficulty with the ISM absorption 
may also influence the accuracy of the bolometric luminosities 
of the stars under investigation. 

The numerical parameters of the S06 scaling relations are 
generally similar to ours: e.g. for their full sample of O- 
type stars and in the total energy band 0.4-lO.OkeV their 
log(Lx/Lboi) = -6.865 ± 0.186 (cf. our value for the 0.4- 
lO.OkeVband: -6.58 ± 0.79). 

The biggest difference between Sana et al. and our values of 
the scaling parameter for O-type stars is in the medium band. It is 
explained by the fact that in our fits we included all O-type stars 
while S06 provide the best fit value for a sub-sample of O-type 
stars excluding two deviating binary systems. If we exclude two 
anomalous stars Tr 1421 and CPD— 59° 2629 from the medium 
band fit, the scaling parameter changes from —7.18 to —7.32, a 
value closer to the one of Sana et al. (—7.58 ± 0.2). Our results 
on the Lx - iboi relation can be summarised as follows: 



1. The canonical relation Lx ^ 10^^ x Lboi is mainly de- 
fined by relatively soft X-ray fluxes below 2.5 keV. In the 
hard band (E > 2.5 keV), the measured X-ray flux is proba- 
bly poorly constrained, except for the brightest sources. As 
a consequence, the scatter of the hard X-ray luminosities is 
large and no definitive fit possible. 

2. S06, based on uniform X-ray data and homogeneous obser- 
vations, state that the observed dependence between intrinsic 
Lx and Lboi is much tighter than was thought before, with 
the residual dispersion of the Lx - iboi relation around the 



The X-ray bands of S06 are almost identical to ours. 
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least-square scaling law being less or equal to some 20%, if 
one excludes binary stars and stars later than 09. Our current 
results are unable to confirm this statement due to a larger 
scatter in the data owing to the patchy ISM absorption in the 
Carina region. 
3. On average, the observed X-ray properties of O binary stars 
are not very different from those of single stars (cf. our 
Fig.[T4]and Fig. 16 from S06). We shall discuss possible rea- 
sons for this behaviour below. 



7.3. On the origin of tlie Lx -iboi scaling for single stars 

It has been pointed out many times in the literature that the scal- 
ing relation between Lx and Lboi could be an indirect effect. 
Indeed, the common paradigm for the generation of X-rays in 
early-type stars is a wind shock model (see e.g. Lucy & Solomon 
[19701 Owocki & Cohen [19991 Dessart & Owocki[2005Tl. If X- 
ray emission is formed throughout the wind then Lx should be 
proportional to some combination of wind parameters like the 
mass loss rate M and the terminal wind velocity Woo. Indeed, 
Owocki & Cohen (,1999 J from the exospheric approximation es- 
tablished that the X-ray luminosity should scale with the mean 
wind density M/uoo in the form Lx ^ (Af/uoo)^ for optically 
thin winds and Lx ^ (Af /woo)'-^^'*-' for optically thick winds, 
with s being the index of the radial power law dependence of 
the X-ray filling factor f ^ r^ . However, Sana et al. ( I2006I I 
state that these considerations contrast with the observed tight 
dependence between Lx and Lboi- 

A new idea about the origin of X-rays from single O-type 
stars was recently brought forward by Pollock (2007). It is based 
on the fact that, according to recent XMM-Newton high res- 
olution observations of (^Orionis, all emission lines have the 
same velocity profile which contrast with the paradigm of wind 
shocks, in which line profiles should be dependent on the global 
kinematic structure of the wind. The principal difference be- 
tween the two models is that in the traditional wind shock model 
the energy exchange between post-shock ions and electrons oc- 
curs on a short time scale while Pollock estimates the ion-ion 
collisional mean-free path and find that at least in the case of 
(^ Orionis it can be quite large, which leads to very slow energy 
exchange between hot ions and cold electrons. He further sug- 
gest that the shocks are collisionless and that the ionization of 
the post-shock gas is not caused by electrons as in the traditional 
picture but rather by protons. The line profiles then are defined 
not by macroscopic motion of ionized gas but simply by the line- 
of-sight component of the thermalized motion of ions in the im- 
mediate post-shock gas. This idea has not yet been developed up 
to a stage when numerical modeling could be done and a theo- 
retical connection between Lx and Lboi (if any) established. 

From the observational point of view, we agree with S06 that 
whatever the theoretical explanation behind the Lx - iboi rela- 
tion, it is a firmly established observational fact which should 
be explained by a future theory. Not only is the proportionality 
itself now well established but also some of its more detailed 
characteristics (see items 1-3 above). 



7.4. Binary versus single stars 

Item (3) in the above list requires special consideration. It was 
already mentioned in the introduction that, in a binary system 
consisting of two early-type stars, the supersonic winds of the 
components will inevitably collide. The key feature of the wind- 
wind collision, making it different from localized shocks formed 



by internal line-driven instabiUty of supersonic winds is that the 
relative velocity of the two flows is very high, of the order of 
thousands kms^^, at least along the line connecting the centers 
of the components. For comparison, the velocity differences in 
localized shocks within a wind are typically of the order of a 
few hundred kms^^. As a result. X-ray emission from binary 
early-type stars should be more luminous and harder than in 
single stars. This was understood long ago (Prilutskii & Usov 
119761 Cherepashchuk [T976l l. One of the first theoretical predic- 
tions about the hardness and X-ray luminosity of binary early- 
type stars was made by Usov ( |1992t and later on was refined 
and quantified in numerical gas-dynamical simulations (see, e.g. 
Stevens & Pittard 1999) . 

Over the years, it became clear that the original somewhat 
naive expectations of more luminous and hard X-ray radiation 
from colliding wind binaries had to be significantly adjusted. A 
large number of factors play a role in the X-ray output of a bi- 
nary, such as the orbital period (hence the distance between the 
components, the pre-shock velocity and the density in the colli- 
sion zone), the wind momentum ratio, the orbital inclination, the 
eccentricity, etc. These factors influence both the rate of intrin- 
sic X-ray production in the collision zone and the absorption of 
the X-rays within the winds of the components. For additional 
discussion, see a recent paper by Linder et al. (120061) . 

It is nevertheless still puzzling that so few early-type bina- 
ries are significantly more luminous and harder than the total 
sample of O-type stars (e.g. in our Fig.[T4land in Fig. 16 of S06, 
also see Raassen et al., 2003| who present a collection of X-ray 
properties of single and binary early-type stars). Two intrigu- 
ing questions already mentioned in the introduction are (i) why 
some known colliding wind binaries have soft X-ray spectra and 
(ii) why some stars having hard X-ray spectra (as demonstrated 
e.g. by the presence of the Fe XXV-XXVI Unes at --6.7 keV) 
lack any evidence of binarity? 

The answer to the first question may be two-fold. First, it 
may be partly the problem of the detection of the hard X-ray 
spectral tail in relatively faint colliding wind binaries. On the 
other hand, theoretical studies of wind-wind colUsion mecha- 
nisms (Stevens & Pollock [T994l Owo cki & Gavlev [T995l Walder 
& Folini 2003 Antokhin et al. |2004| l suggest that either the so- 
called radiative braking or radiative inhibition (mutual radiative 
influence of the binary components on the wind of their vis-d- 
vis) may significantly slow down the winds of the binary compo- 
nents just before the collision, thus reducing the X-ray luminos- 
ity of the system and the hardness of its X-ray spectrum (which 
is proportional to the square of the wind velocity entering the 
collision zone). To verify this possibility, we shall apply the re- 
cent models accounting for this phenomenon, in our forthcoming 
analysis of the binary HD 93205 which has one of the best X-ray 
coverages of the orbital cycle of any early-type binary. 

There is probably no single answer to the second ques- 
tion, either. In some cases (most notably WR25 and MJ596) 
it was recently found (see references above) that these hard X- 
ray sources are indeed binary systems that escaped detection be- 
fore. These recent advances raise interest in a similar case of 
a bright 03 V star HD 93250 (showing clear indication of the 
Fe XXV-XXVI lines). The star does not display any radial ve- 
locity or photometric variability and as such was considered as 
single. However, it was reported to display non-thermal radio 
emission (Leitherer et al. 1 19951 ) and this feature is nowadays of- 
ten believed to be evidence for a colliding wind system also in 
Oh-O binaries. On the other hand, the 08.5 V star CPD -59°2629 
which is also an unusually bright and hard X-ray source for its 
spectral type might be an example of a star with a magnetically 
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confined wind (ud-Doula et al. |2006i A low-mass companion 
with high activity could also be considered. 

8. Conclusions 

High quality and large amounts of XMM-Newton data obtained 
for the Carina Nebula, allowed us to detect 235 discrete sources 
in the field of view. Several of these sources are probably pre- 
main sequence stars with a characteristic short-term variabil- 
ity, while seven sources are recognized as possible background 
AGNs. 

We concentrated in our study on the properties of early-type 
stars in the region. Spectral analyses of twenty three sources 
of type OB and WR 25 were performed. We derived the spec- 
tral parameters of the sources and their fluxes in three energy 
bands. Estimating a value of interstellar absorption towards ev- 
ery source and assuming a distance of 2.5 kpc to the Nebula, 
we derived the X-ray luminosities of these stars. Comparison of 
these luminosities with bolometric luminosities has confirmed 
the known scaling relation Lx ~ 10^^ x Lboi for early type 
stars and allowed us to make a few more detailed statements 
about the characteristics of this relation. We pointed out that, on 
average, the observed X-ray properties of binary and single early 
type stars are similar, and give several possible explanations for 
this fact. 

Further work will be aimed at the detailed analysis of the 
data of particular sources. 
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